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ABSTRACT 



Context. This paper is part of our ongoing study of star formation in the (far-) outer Galaxy. 

Aims. Our goal in this paper is to study stars and molecular gas in the direction of IRAS06145+1455 (WB89-789). The kinematic 
distance of the associated molecular cloud is 11.9 kpc. With a galactocentric distance of ~ 20.2 kpc, this object is at the edge of the 
(molecular) disk of the Galaxy. 

Methods. We use near-IR (J, H, K), molecular line-, and dust continuum observations. 

Results. The near-IR data show the presence of an (embedded) cluster of about 60 stars, with a radius ~ 1.3 pc and an average stellar 
surface density ~ 12 pc~ 2 . We find at least 14 stars with NIR-excess, 3 of which are possibly Class I objects. The cluster is embedded 
in a ~ 1000 M G molecular/dust core, from which a molecular outflow originates. The temperature of most of the outflowing gas is 
*L40 K, and the total mass of the swept-up material is *! 10 M Q . Near the center of the flow, indications of much higher temperatures 
are found, probably due to shocks. A spectrum taken of one of the probable cluster members shows a tentative likeness to that of a 
K3 Ill-star (with an age of at least 20 Myr). If correct, this would confirm the kinematic distance. 

Conclusions. This cluster is the furthest one from the Galactic center yet detected. The combination of old and recent activity implies 
that star formation has been going on for at least 20 Myr, which is difficult to understand considering the location of this object, where 
external triggers are either absent or weak, compared to the inner Galaxy. This suggests that once star formation is occurring, later 
generations of stars may form through the effect of the first generation of stars on the (remnants of) the original molecular cloud. 

Key words. Stars: formation - Stars: pre-main sequence - ISM: clouds - ISM: individual objects: WB 89-789 (IRAS06145+1455) 



1. Introduction 

The Hi in our Galaxy extends out to galactocentric distances R 
of at least 24 -25 kp c (e.g., Wouterloot et al. I1990I McClure- 
Griffiths et al. 2004). From high-sensitivity 21 -cm observations, 
Knapp et al. ( 1978) found Hi out to R ~ 50 kpc (for a flat ro- 
tation curve and Ro=8.5 kpc). Molecular material and associ- 
ated sites of star formation do not appear at such large distances, 
however. Star formation in the Galaxy has been observed to oc- 
cur out to R ss 20 kpc (e.g., Fich & Blitz [1984 1; Wouterloot et 



al. 1988; Kobayashi & Tokunaga 2000; Santos et al. 2000; Snell 
et al. 120021) . Wouterloot & Brand (1989) detected many molec- 
ular clouds (hence star formation sites) towards IRAS sources 
in the outer regions of the disk. Thus both star formation and 
star formation reservoirs (i.e., molecular clouds) are present at 
large distances from the Galactic center, but not as far out as the 
atomic hydrogen. This suggests that beyond about R « 20 kpc, 
conditions are not favorable to ttansforming Hi into H2 and then 
into stars. The cause of this is not known (in general the for- 
mation of H2 clouds from Hi is not yet a settled manner, but it 
must have to do with the physical environment presented by the 
interstellar medium). 

In the far-outer Galaxy (R > 16 kpc; hereafter FOG) the 
physical environment differs from that in the inner Galaxy (for 
a summary see Brand & Wouterloot |19951 l, the effects of which 
may influence the formation of molecular clouds and the star 
formation process within them. Also, both the volume density 
of the molecular (and atomic) gas and the strength of spiral den- 
sity waves are much reduced in the outer Galaxy. Even if star 



formation is independent of cloud formation and is to be initi- 
ated or enhanced by cloud-cloud collisions or triggered by su- 
pernovae or density waves, then star formation activity in the 
outer Galaxy would be expected to be significantly lower than 
in the inner Galaxy. To study the process of star formation and 
its end products in the outer parts of the galactic molecular disk, 
we observed a selection of FOG clouds in the NIR to directly 
detect the embedded stellar population. FOG clouds are at large 
enough R for the physical conditions of the ISM to be different, 
yet those in the 2 nd and 3 rd quadrants are close enough for the 
star-forming cores to be resolved. 
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Here we present the most distant object studied with 
these observations: WB89-789 (IRAS06145+1455; a 20 oo = 
06 /, 17" ! 2452, (5 20 oo=+ 14°54'42") at (l,b = 195.°82, -0.°57), 
which is located at R « 20.2 kpc, d « 11.9 kpc (Brand & 
Wouterloot l 19941 hereafter BW94). The luminosity of the IRAS 
source (Lfi r ~ 1.9 x 10 4 L ) is consistent with a single star of type 
B0.5V; it is embedded in a small (equivalent radius r e m 5.7 pc, 
and M ~ 5.6 x 10 3 M [BW94]) cloud. A n H 2 Q maser was 
detected towards this object (Wouterloot et al.[T993), indicating 
that this is an active star-forming region. No radio continuum 
emission was found within 4' of the location of the IRAS source 
in our VLA A-array observations at 3.6 cm (Nov. 1993; unpub- 
lished). With a Act detection limit of ~ 3 mJy we would have 
been able to detect an Hn region ionized by a B0.5V star or ear- 
lier. 
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2. Observations and data reduction 

2.1. Near-infrared 

The data presented here were obtained on February 15, 1995, 
with the ESO 2.2-m telescope at La Silla (Chile). Images in J, H, 
and K-bands were taken with the IRAC-2 camera and objective 
C, which resulted in a scale of 0."49/pixel. The detector was a 
256 x 256 pixel 2 NICMOS-III. The total instantaneous field-of- 
view is therefore ~ 2.1 x 2.1 arcmin 2 . During the observations 
the seeing was 0."9. 

From a quick K-band exposure (20 x 3 seconds) centered on 
the WB 89-789 IRAS point source position, a group of stars was 
seen to be small enough in extent to allow an observing proce- 
dure where it was placed at the center of each quadrant of the 
CCD, as well as near the center of the CCD itself. Integrations 
in the J-, H-, and K-bands were then made for 60 x 2 seconds, 
resulting in 5 separate images per band. After correction with a 
bad-pixel mask, each image was reduced separately by subtract- 
ing a starless sky image derived by median averaging from the 
other 4 exposures, and dividing by the flat field (see below). The 
reduced images were then averaged to get the final frame, which 
has an effective integration time of 10 minutes. Immediately af- 
ter these "ON-source" observations a comparison (OFF-) field 
was observed, at the same galactic latitude, but shifted towards 
smaller longitude by 0.°2 (* 42 pc at the distance of WB89-789). 
To construct flat fields, dome exposures in all three bands were 
taken at the start and at the end of each night. 

We observed 5 standard stars during the night, at different 
airmass. Each standard star was observed at 5 different positions 
on the CCD (with the same observing pattern as used for our 
object) in all 3 bands. Each observation was reduced individually 
and the derived photometric zero points (ZPs) were averaged, 
after which the ZP at the airmass of the WB89-789 observation 
was obtained through linear interpolation. 

The standard star observations revealed a non-uniform CCD, 
with some quadrants "hotter" (i.e., registering more electrons for 
the same number of infalling photons) than others. This non- 
uniformity will increase the random error in the photometry of 
the science objects. In part this is accounted for by the average 
ZP, as the standards were observed in the same way as the sci- 
ence object, hence the average ZP partly takes into account the 
differences between the quadrants. 

The mosaiced images in J, H, and K were aligned and trun- 
cated to show the same region of the sky. Stars were identified 
in each frame with the IRAF task DAOFIND, and photometry 
was performed through DAOPHOT, using a point-spread func- 
tion determined from 5 relatively bright and isolated stars. Fitted 
stars were subtracted from the images, after which DAOFIND 
was run again on the residual frame. This procedure was re- 
peated until no more stars were found (usually after 2 or 3 it- 
erations). Derived magnitudes were then corrected to take into 
account the difference between the point spread fitting radius (3 
pixels) and the aperture radius used (30 pixels) on the standard 
stars to derive the ZP; this correction amounts to -0!"34, -0!"27, 
and -0I"24 for J, H, and K, respectively. 

2.2. Millimeter-lines 
JCMT 

On January 24, 2003, we used the 15-m James Clerk Maxwell 
Telescope (JCMT) on Mauna Kea (Hawaii) to map the emission 
of 12 CO(2-l) towards WB89-789 in the raster (on-the-fly) ob- 
serving mode. The JCMT beam size at 230 GHz is 22"; obser- 



vations were made on a 10" raster, using an autocorrelator spec- 
trometer at a velocity resolution of 0.2 kms -1 . An off-position 
was used at offset 600" in right ascension from the IRAS source 
position. The typical rms noise level in the spectra is 0.2 K (T^). 

On August 7 and 8, 2003, the JCMT was used to simulta- 
neously observe 13 CO(2-l) and C ls O(2-l) emission in a 80"x 
80" region on a 10"-raster towards WB89-789 using frequency- 
switching over 8.2 MHz. The typical rms noise level in these 
spectra is 0.13 K (TV). A comparison with separate (i.e., non- 
simultaneous) 13 CO(2-l) observations of line calibrators shows 
that our 13 CO(2-l) intensities must be corrected by multiplica- 
tion with a factor of 0.88. 

On January 5, 2006, we used the JCMT to observe the 
12 CO(3-2) transition towards WB89-789, in essentially the same 
region covered by the 12 CO(2-l) observations. Observations 
were made on a grid with 6" spacing; the beam size at 
345.796 GHz is 14". The velocity resolution is 0.136 kms" 1 , 
and typical rms noise level is 0.6 K (TV). Pointing was found to 
be accurate within 1 ."8. The line-data were reduced and analyzed 
with the programs CLASS and GRAPHIC, which are part of the 
GAG-software package developed by the Obs. de Grenoble and 
IRAM-Grenoble. 

IRAM 30-m 

WB89-789 was observed in CS J=2-l, 3-2, and 5-4 with the 
IRAM 30-m telescope at Pico Veleta (Granada, Spain) in the pe- 
riod July 24 - 26, 1991. Maps were made of thirteen positions, 
with 15" spacing. The region mapped was sufficient to cover all 
CS emission. The telescope beamwidth at these frequencies (98, 
147, and 245 GHz) is respectively 26", 16", and 11". All intensi- 
ties are on a T* K scale. For CS(2-1) and (3-2) we used a 2x128- 
channel filterbank of 100 kHz per channel, resulting in a reso- 
lution of 0.31 and 0.20 kms~', respectively; CS(5-4) was ob- 
served with a resolution of 584 kHz (0.71 kms -1 ). Observations 
were made by position switching against an off-position 30' W; 
the typical rms noise level in the spectra at the central position 
is 0.05 K. All spectra were reduced with the Grenoble CLASS 
software. 

2.3. Millimeter-continuum 
JCMT 

On April 10, 2002, we simultaneously observed the 450 fim and 
850 fim emission towards WB89-789 with the Submillimetre 
Common-User Bolometer Array (SCUBA; Holland et al. I1999D 
at the JCMT. A standard 64-point jiggle map was made, which 
covers an area 2'3 in diameter; hence the observations are 
limited to the region immediately around the IRAS source. 
Calibration was performed by using fits to the time-dependence 
of skydips and the nearby Caltech Sub-millimeter Observatory 
(CSO) 225 GHz optical depth data, and average ratios of 
r(450 fim) and r(850 /urn) to r(225 GHz). The final transi- 
tion from instrumental parameters to Jy/beam was derived from 
similar maps of Mars. The half-power beam size at 450 fim 
and 850 fim are 9" and 15", respectively. The rms values 
in the WB89-789 maps are 450 mJy/beam (at 450 fim) and 
20 mJy/beam (at 850 fim), respectively. The data were reduced 
with the program SURF (Jenness & Lightfoot 1998). We found 
that away from the continuum source, the background emission 
level at both wavelengths was not around zero, but somewhat 
higher. We have corrected this by subtracting 1 Jy and 0.05 Jy 
from all pixel values in the 450 fim and 850 fim maps, respec- 
tively. 
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SEST 

Between August 17 and 22, 2003, the region around WB89- 
789 was observed in the 1.2-mm continuum with the 37-channel 
bolometer array SIMBA (SEST IMaging Bolometer Array) at 
the 15-m Swedish-ESO Submillimetre Telescope (SEST) at 
ESO-La Silla, Chile. A region of size 600" x 600" (azimuth 
x elevation) was scanned at a rate of 80"/s; the total integration 
time per map was about 12 minutes. Atmospheric opacity was 
determined from skydips, which were taken every 2 hours; val- 
ues at the zenith ranged between 0.15 and 0.45. The data were 
transformed from counts/beam to mJy/beam by making similar 
maps of Uranus once per day. The average conversion factor dur- 
ing this observing period is 58+6 mJy/count. 
In total, 8 maps were obtained. Each map was reduced and cal- 
ibrated individually, after which they were averaged. The rms 
in the final image is 15 mJy/beam. Pointing of the SEST was 
checked by observing a strong continuum source every 2 hours. 
The half -power beam width of the SEST at 1.2 mm is 23". Data 
were reduced with the program MOPSI, written by Robert Zylka 
(IPvAM, Grenoble), following the instructions from the SIMBA 
Observers Handbook (Version 1.9, 9 Feb. 2003). 

2.4. Spectroscopy 

On December 7, 2004, we used the DOLORES spectrograph at 
the Telescopio Nazionale Galileo (TNG) at La Palma (in service 
observing-mode) to take a spectrum of one of the stars of the 
cluster. Six separate integrations of 2220 sec each were taken us- 
ing the high-resolution spectrograph and HR-B Grism nr. 5 (dis- 
persion 0.875 A per pixel), which resulted in a wavelength cov- 
erage of about 3200A-4900A, and (with a 1" slit) in a spectral 
resolution of ~3.2A. Bias frames were taken on the same night, 
and flat fields and He and Ar calibration-lamp spectra were taken 
the night before and after. As a flux calibrator, Hiltner600 was 
observed on the same night. 

Data reduction was performed in IRAF, using standard pack- 
ages. Bias frames were subtracted from all science and flat 
field frames. An average normalized flat field was constructed, 
which was corrected for bad pixels and then divided into the 
star frames. Wavelength calibration was performed with the Ar- 
lamp spectra, which in the wavelength range in question only 
has 8 useful (i.e., non-blended) lines, none below 3780A, and is 
therefore not very accurate. 



3. Results 

3.1. Stars 

3.1.1. Photometry, clustering, and pre-main sequence stars 

In Fig.QJ we show the K-band mosaic, while a false-color (JHK) 
image of the area around WB89-789 (the region inside the box 
in Fig.QJ) is shown in Fig.QJ). Just below the center of the image 
we see a group of stars embedded in nebulosity; the IRAS source 
is located in the northern half of the nebulosity. In this 1 .'7 x 1 .'7 
area DAOFIND detects 173, 179, and 161 stars in the J-, H-, and 
K-bands, respectively. The distributions of these stars in bins of 
1 magnitude peak at raj ~ 20 (faintest star found 19.96±0.24), 
m H ~ 18 (19.59+0.54), and m K * 17 (18.52+0.50), respectively. 
Signal-to-noise ratios are reached of 5-10, i.e., photometric ac- 
curacy of 0!"l-0!"2, at /«k = 16.5 - 17.5, which corresponds to 
main sequence [ms] spectral types A0-5 (2-3 M ) for a distance 
of 10 kpc and visual extinction Ay = 10 mags. 
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Fig. 2. Comparison between the star densities in the WB89-789 
field (drawn) and in a field centered at a nearby off-position 
(dashed). Stars were counted on the K-images in concentric rings 
of equal area. In the WB89-789 field the rings were centered on 
the nebulosity (Fig. QJi). An excess of stars in the WB89-789 
field is clearly visible. 



To confirm that WB89-789 is associated with a clustering of 
stars, in Fig. [2] we compare the density of stars in the WB89- 
789 field with that at a nearby off-position (see Sect. [2]). Stars 
were counted in concentric rings of equal area (n x 15(pixels) 2 ). 
For the WB89-789 field the rings were centered on the nebulous 
region visible in Fig. [3] and the counting was performed in the 
boxed-in area of the K-frame shown in Fig. [T^. There is a clear 
excess of stars compared to the background out to about 22" 
from the center of the nebulosity; this corresponds to a radius 
of 1.3 pc, in which the average stellar surface density is 12 pc~ 2 
(compared to <3 pc~ 2 for the OFF-field). There is a particularly 
strong concentration of stars in the inner 15 pixels (~ 7."5); in 
this region, with r ~ 0.4 pc, the stellar surface density is 33 pc -2 . 

The central 0.'9 x l.'l part of the region shown in Fig. [TJ? is 
where all 5 individual ON-source images overlap; a (K-frame) 
contour plot of this region is shown in Fig [3] Sixty-eight stars 
were detected in all three (J, H, K) bands. The photometric data 
are collected in Table[4] The (J-H), (H-K)-diagram of these stars 
is shown in Fig. [4] The dotted lines define the reddening band for 
normal stellar photospheres (Rieke & Lebofsky 1985). Objects 
outside and to the right of this band have intrinsic NIR excess, 
and for a number of them their location in this diagram can be 
explained by the presence of circumstellar disks, with central 
holes of various dimensions (see Lada & Adams [19921 ); these 
are likely pre-ms stars. 

There are 4 stars outside the normal reddening region on the 
left (nrs. 87, 53, 51, and 20). All of these except nr. 87 are consis- 
tent with a location inside the normal region, considering error 
bars. On the other hand, stars with effective temperatures be- 
tween 2500 K and 4000 K are predicted to have an excess of 
emission in the H-band, which is due to the formation of H2 and 
absorption by H (Gingerich & Kumar 1964); this will shift stars 
like nr. 87 to the left of the normal reddening region in Fig.|4](see 
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Fig. 1. a (left). K-band image of the region around WB89-789 (IRAS06145+1455). The area of sky visible here is ~ 3 x 3 arcmin 2 . 
The region outlined by the black box has a size of ~ 1 .7 x 1 .7 arcmin 2 , and is shown on the right. North is up, East is left, b (right). 
False-color image (J=blue, H=green, K=red) of a ~ 1 .7 x 1.7 arcmin 2 region around WB89-789 



also Chini et al. 119921 . There are 22 stars outside the normal re- 
gion on the right. Eight (nrs. 7, 9, 13, 40, 63, 81, 98, 99) have 
anomalous colors, in that they are also below the unreddened 
TTau locus; these are indicated with circles in Fig. [6] Three of 
these (nrs. 13, 40, 98) could be brought into normal region con- 
sidering their error bars. The other 5 have 'anomalous' colors, 
which could be due to unresolved binaries (Lada et al. 2000). 
That leaves at least 14 (22-8) stars with 'true' NIR excess (i.e., 
21% of the total of 68); these are indicated with a star in Fig. [6] 
Eleven are in the Class Il-source region (between the right-hand 
normal reddening line and the right-most TTau reddening line). 
Their SEDs are dominated by disk emission and their location 
can be explained by variation in size of central hole and incli- 
nation (see Lada & Adams 1992). Three (nrs. 15, 54, 83) are in 
the region of Class I sources. These are the most embedded and 
youngest objects. Their SED is dominated by emission from an 
envelope of gas and dust. 



In Fig. |5j3 we show all 68 stars for which we have J, H, and K, as 
well as the 15 stars for which only H and K are available. Again 
ignoring stars with NIR-excess (identified with a large open cir- 
cle), as well as the stars for which we only have H and K (for 
which we cannot say whether they have NIR-excess) we see that 
there are stars with an additional extinction of up to ~ 30 mag- 
nitudes, although the bulk has Ay ~ 7.5. This is about the same 
for stars inside and outside the above-defined cluster radius. As 
we mentioned in the Introduction, the edge of the molecular disk 
of the Galaxy lies at about R w 20 kpc, thus it is unlikely that in 
our images of WB89-789 there are any stars more distant than 
the cluster, and we conclude that most stars shown are part of the 
cluster. Potentially the earliest-type star in the cluster that does 
not have NIR excess is star 33 (see Fig[3j: If one were to simply 
shift it back to the the main sequence, it would be of spectral 
type B0, consistent with the spectral type derived from the 
(and assuming it is an ms-star). 



3.1.2. Extinction 



The amount of visual extinction towards WB89-789 is de- 
rived by means of a K versus (H-K) color-magnitude diagram 
(Fig. |5^), plotting only those stars that are within 7 ."5 from the 
center of the nebulosity and thus presumably part of the cluster. 
Here we can also use those (7) stars that were detected only in 
H and K, but not in J. In addition to the stars, we have drawn 
the location of the unreddened main sequence for a distance of 
11.9 kpc. Ignoring the stars that have been identified as having 
a NIR-excess in Fig. [4] because they do not follow the standard 
interstellar reddening law, we shift the main sequence along the 
direction of the interstellar reddening (indicated by the arrow in 
Fig. [5] A K = 0.112 x A v and E(H-K) = 0.061xA v ; Rieke & 
Lebofsky 1985) until it reaches the locations of the cluster stars. 
Thus we derive a maximum amount of foreground visual ex- 
tinction of Ay ~ 6.25 magnitudes (maximum, because even the 
least-reddened stars can have some amount of internal extinc- 
tion, i.e., due to the dust in the molecular cloud in which they 
are embedded). This may not seem much for such a supposedly 
distant object, but it is consistent with the findings of, e.g., Fich 
( 1 1984| l and Amores & Lepine (2005), who concluded that the 
extinction is fairly low (~ 3 mags.) in much of the outer Galaxy, 
and that there does not appear to be a progressively increasing 
extinction beyond a few kpc from the Sun. 



Figure|6]shows a section of the K-band image of the WB89- 
789 area, with the contours of integrated C 18 0(2-l) emission 
(see Fig. [7};) superimposed. The 14 stars with NIR-excess are 
shown in Fig. |4j), and the 8 stars with anomalous colors are 
marked in Fig. [6] Rather than being concentrated in the center 
of the cluster (which we take to be the nebulous region seen in 
Fig. [3}, the NIR-excess stars are distributed in a ring around the 
peak of the C ls O(2-l) emission. At the location of the molecular 
peak, no stars are seen in this K-image, which may be because 
the extinction there is too high. The average column density of 
H 2 , N(H 2 ), derived from the C ls O(2-l) data (see Sect. 3. 3.1) is 
7.8xl0 21 cm 2 . Because N(H 2 )=10 21 xA v (Bohlin et al. [19781 
and assuming the gas-to-dust ratio in the FOG is the same as lo- 
cally) this implies a visual extinction Ay of about 8 magnitudes 
due to the whole of the molecular clump (front to back). At the 
peak of the C 18 0(2-l) emission N(H 2 )« 1.0 x 10 22 cirT 2 , and 
thus Ay « 10 magnitudes. If the cluster stars are embedded in 
the molecular gas, as they are likely to be, then the general (i.e., 
not counting that which is due to the circumstellar material) in- 
ternal extinction would therefore be on average ~5 magnitudes 
(and in any case <10 magnitudes), consistent with the number 
derived above from the color-magnitude diagram. 
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Fig. 3. a. Contour plot (K-frame) of that part of the final mosaic where all 5 ON-source images overlap. Sixty-eight stars in this 
area have been detected in all three bands (J, H, and K). Several stars are labeled, to identify the field in panel b. Among them is 
star 33, of which a spectrum has been taken (Sect. 2.4). North is up, East is left. The (0,0)-offset indicates the nominal position of 
the IRAS source. 

b. The cluster field from panel a (inside the box) shown in a larger context on the POSS ITR image. 
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Fig. 4. a. NIR color-color diagram for the 68 stars in Fig. [3^ that have reliable detections in all three (J, H, and K) bands. The thick 
and thin drawn curves are the unreddened main sequence and giant branch, respectively, from Bessell & Brett (119881 1; the dashed 
line indicates the loci of the classical T Tauri stars (Meyer et al. 11997k the dotted lines indicate the direction of normal interstellar 
reddening (Rieke & Lebofsky 1985). Crosses on these lines mark increments of 5 magnitudes of visual extinction from the points 
where they intersect the main sequence curve. Objects outside and to the right of the reddening band have intrinsic NIR excess; to 
avoid confusion, error bars are shown only for objects that are well outside the normal reddening band. The filled triangle identifies 
star 33 (see text and Fig. [3). b. Like a, but with the stars identified by their ID-number. 



3.2. Molecular cloud 
3.2.1. Cloud parameters 

Whole cloud 

BW94 mapped the molecular cloud associated with WB89- 
789 in CO(l-O) with the SEST (beam size 45"), on a 40" raster. 
We have obtained new, better sampled CO(2-l) and (3-2) maps 
with the JCMT; the latter is shown in Fig.EJi. From the CO(2-l) 
data, following the procedures outlined in BW94, we derive a 
mass M co ~ 4.5 x 10 3 M (from N(H 2 )=X j T^(CO)dv, with 



X = 1.9 X 10 2() crrT^Kkms -1 ) -1 and using r mb = / r*/^ fss dv; 
?7f ss = 0.80 is the forward efficiency). The virial mass M v „ » 
5.1 x 10 3 M (M vir = 126r e (AV) 2 for an r~ 2 density distri- 
bution; MacLaren et al. 119881 1; AV is the FWHM line width, 
and the equivalent, beam-corrected radius r e « 5.3 pc. These 
numbers agree very well with those derived by BW94 from 
the CO(l-O) map. We also note that the virial mass and the 
mass derived from the empirical method, using X, are in ex- 
cellent agreement. The value of X quoted above is the inner 
Galaxy value; Brand & Wouterloot (1995) have shown that in 
the outer Galaxy X may be somewhat larger, but lies within 
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Fig. 5. K versus (H-K) color-magnitude diagram for WB89-789. a. Only stars within 15 pixels (~7."5) from the center of the nebula 
(Fig.QJ). Filled circles are stars detected in all three bands (J, H, and K); those with a NIR-excess have been additionally marked 
with an open circle. Open triangles: stars detected in H and K only. The drawn line marked 06-M5 is the unreddened main sequence 
for a distance of 1 1 .9 kpc; data from Koornneef (1983] 06 - B7) and Bessell & Brett ( 1988 ; B8 - M5). K-magnitudes were derived 
from My (Schmidt-Kaler 1982) and (V-K)q. The dashed line is the main sequence for a distance of 1 1.9 kpc and a visual extinction 
Ay = 6.25 mags (Rieke & Lebofsky 1985 ). The drawn line marked G0-M5 is the unreddened giant branch (data from Bessell & 
Brett 1988 ) for a distance of 1 1 .9 kpc. Its location for Ay = 6.25 mags is shown as a dashed line parallel to the unreddened sequence. 
Star 33 (see text and Fig.0 is indicated, b. Like a, but for all stars. 
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Fig. 6. K-band image of the WB89-789 region, with contours of 
the integrated C 18 0(2-l) emission superimposed (see Fig.fTJ; for 
details). The IRAS position is at (0,0). The stars with NIR-excess 
are labeled, and identified with asterisks; those with anomalous 
colors (see text) are marked with circles. 



30 -45% of the inner Galaxy value. Hence at the edge of the 



Galaxy, X may be 2.8 x 10 20 cm 2 (Kkms ') which results in 
M C o ~ 4.5-6.6 x 10 3 M . These results are collected in TableQ] 

Cloud core: CO 

In Fig. |7j3 and c we show a plot of the integrated emission of 
I3 CO(2-l) and C I8 0(2-1), respectively. The emission in these 
two isotopomers was mapped only in the core of the cloud, but as 
can be seen from the maps, at least for C I8 0(2-1), all the emis- 
sion is contained within the mapped region. While the I3 CO and 
I2 CO(2-l) emission peak at the IRAS position (offset 0",0"), the 
C ls O(2-l)-peakis offset by about 10" to the East; this is not due 
to pointing uncertainties because I3 CO and C 18 were observed 
simultaneously. The dust continuum emission (see Fig. ITOb also 
peaks very close to the IRAS position; hence the offset detected 
in the C I8 emission may be due to that molecule's depletion 
near the core in which the IRAS source is embedded (cf. Tafalla 
etal. 12004] ). The 12 CO emission, being more optically thick, only 
traces the outer layers of the core; also, 13 CO is more optically 
thick than C 18 (see later) and traces the more outlying layers. 
The C ls O(2-l) core is elongated in a SW-direction, similar to 
the 12 CO core in panel a. The core-radius derived from the I3 CO 
observations is identical to that, derived from 12 CO, while the 
C ls O core is ~40% smaller. The masses Mqo and M v ; r deter- 
mined from 12 CO in the same region as mapped in 13 CO and 
C ls O are reported in Cols. 4 and 6, respectively, of TableQ] un- 
der " 12 CO, core". 

Column densities N can be derived from r m b( 13 CO) and 
7 , mb(C 18 0), assuming the emission is optically thin (see, e.g., 
Rohlfs & Wilson [19961 Brand & Wouterloot H^98l Eq. 1). The 
excitation temperature T ex is derived from the optically thick 
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A& (arcsee) Act (arcsec) V, sr (km/s) 

Fig. 7. JCMT-data. a. Map of the integrated intensity of CO(3-2) associated with WB89-789 ( f T* A dv; the integral is over the bulk 
of the emission, 30 < Vi sr < 38 kms -1 , and excludes the line wings). Contour levels 2.5(2.5)10(5)55 Kkms 1 (lowest(step)highest). 
The filled triangle marks the peak of the SCUBA 850 fim emission (at -l."9, 2."1), and the observed positions are indicated 
with the tiny crosses. The IRAS point source position is at offset (0,0). b. Map of the 13 CO(2-l) emission (area under the fitted 
Gaussians). Contour levels are 1(2)15 Kkms -1 . Triangle and crosses as in a. c. As b, but for the C ls O(2-l) emission. Contour 
levels are 0.3(0.3)2.1 Kkms -1 . d. Line profiles at offset position (10",0"). The strongest emission is the I2 CO(2-l)-line; the profile 
clearly shows self-absorption. The less intense emission, plotted as a histogram, is the 13 CO(2-l) line, and the drawn line is that of 
C 18 Q(2-1), multiplied by 4. e. Like d, but for offset position (0", 10"). 



Table 1. Cloud parameters derived from the molecular line mea- 
surements. 





r\ 

pc 


< 

hp 


AV 
km s 


^co 


M L 
10 3 M G 


M 4 

vir 


12 CO, cloud 


5.3 


1.6 


2.77 5 


4.5-6.6 




5.1 








2.69 56 






4.8 6 


12 CO, core 


2.9 


1.6 


3.34 5 


3.2-4.7 




4.1 








3.22 5 ' 6 






3.8 6 


13 CO, core 


2.7 


1.1 


2.14 




1.0-5.0 


1.6 


C 18 0, core 


1.6 


0.9 


1.72 




0.8-3.2 


0.6 


CS, core 


1.1 


0.7 


2.07 






0.4 



1. r e : Equivalent radius, corrected for beam size, at the 2 - 3<x-level; 
rh p : the radius at half-maximum intensity level, corrected for beam 
size. 

2. ViaN(H 2 )=X / T mb (CO)dv. The indicated range of masses derives 
from using the inner Galaxy value X = 1.9 x 10 20 cirr 2 (Kkm s" 1 ) -1 and 
allowing a possible 45% increase at the edge of the Galaxy (see text). 

3. Assuming LTE-conditions, and T cx =20 K. Smaller values are for 
local abundance ratios for [H 2 ]/[ 13 CO] and [H 2 ]/[C ls O], higher values 
are for abundances extrapolated to the FOG (see text). 

4. M vir =126 r (AV) 2 , for a density distribution oc r~ 2 , and using r e . 

5. Equivalent width, derived from the area under the emission line and 
the peak temperature. 

6. Excluding the wing emission. 



'"CO transition. The CO emission line is not just optically 
thick: the line profiles within 20-30" from the IRAS position 
show self-absorption. Examples are shown in Fig. [7}i, e; we 
shall return to this in Sect. 3.2.2. From the emission at posi- 
tions away from the region with self-absorption, we find that 
r ex ~ 20 K is a representative value for this cloud. At each grid 
position of the cloud-map we derived N( 13 CO) and N(C ls O). 
Column densities of H2, and LTE-masses M\ te were then derived 
through N(H 2 ) ~5x 10 5 N( 13 CO) (Dickman & Clem ens [19831 
and N(H 2 ) ~ 6 x 10 6 N(C 18 O) (Frerking et al. [19821 The M lte 



are presented in Col. 5 of Table [T] From the detection equation, 
and the adopted value of Tqx, we find that on average r(C^O) 
~ 0.05 (peak value 0.11 at offset 10",0"), and t( 13 CO) ~ 0.22 
(peak value 0.64 at offset 0",-10"). It is clear that the C lg O is 
optically thin, whereas the 13 CO is less so. 
The Mite derived from 13 CO and C 18 are in good agreement, but 
they are both smaller than the core-masses derived from 12 CO. 
This is partly due to the fact that in I2 CO one sees the outer layers 
of the cloud, not seen in the other two isotopomers. In the deriva- 
tion of the LTE-masses we have used the locally-determined 
abundance ratios [H 2 ]/[ 13 CO] and [H 2 ]/[C 18 0]. However, be- 
cause of the galactic abundance gradients these values are cer- 
tainly not valid at the edge of the Galaxy. As argued by Brand & 
Wouterloot ( 119951 1, the ratio [H 2 ]/[ 13 CO] in the far-outer Galaxy 
might be about 5 times higher than locally. A similar extrapo- 
lation of the abundance gradients (Wilson & Matteucci 1992; 
Wilson & Rood |1994| l suggests that the ratio [H 2 ]/[C ls O] in 
the far-outer Galaxy might be about 4 times higher than the lo- 
cal value. Taking this into account, M\ te derived from 13 CO and 
C 18 would become 5 x 10 3 M and 3.2 x 10 3 M , respectively, 
which is closer to the range for Mqo (see TableQ}. 

Cloud core: CS 

All CS lines are detected essentially only at the IRAS point 
source position. The lines are Gaussian, showing no signs of 
outflow or other asymmetries. From the CS(2-1) map we find 
a (beam-corrected) radius of the CS-core of 0.7 pc, smaller 
even than the C ls O core. With an average line width AV « 
2.07 kms _I , we derive a virial mass of M wu « 400 M , similar 
to what was found from the C 18 data. We have convolved the 
CS(3-2) and (5^4) spectra at the central position to the CS(2- 
1) beam, assuming Gaussian beam profiles. The resulting spec- 
tra were put on a main-beam brightness temperature scale (T m b), 
and used in an LVG-model to derive density and kinetic tempera- 
ture of the CS-gas. We ran the model with a range of parameters: 
T=10-100 K in steps of 2 K, log(n H2 )=4-6 in steps of 0.1-0.2, 
and abundances lxlO" 11 , 3 x 10" 11 , 5 x 10 -11 , and 1 x 10~ 10 . We 
used a velocity gradient dv/dr ~ 1.5 kms^'pc" 1 . We found the 
smallest^ 2 (~ 0.6), hence the best reproduction of the observed 
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line ratios, for a CS-abundance of 3 xlO~ n , T^ n =s 84 K, and 
log(nH,/cirT 3 ) m 4.9. The fit is not very sensitive to the temper- 
ature, as x 1 1S within a factor of 2 of the minimum value over 
a range of about 30 K. Nevertheless it is clear that the deeper 
layers of the core from which the CS emission originates are 
considerably hotter than the outer parts traced by 12 CO. 

3.2.2. Molecular outflow 

The SEST CO(l-O) spectrum at the IRAS position showed a 
line profile with non-Gaussian wings (BW94). The new, higher- 
resolution CO(2-l) data allow us to better outline the outflow 
(see Fig. [8^), which is completely contained within one SEST 
CO(l-O) beam (~ 45"; Fig.fjp). In Fig. [8^ we show a position- 
velocity diagram at A6=0" . Enhanced emission in the line wings 
is clearly visible at Aa-0" . Also clearly visible is a double-peak 
structure in the 12 CO(2-l) emission. At positions within ~ 20" 
from the center of the map the lines show a significant inden- 
tation, as illustrated in Fig. [7Ji, e. By comparing the position- 
velocity map in Fig. [8^ with the theoretical contour maps shown 
in Figs. 3 and 4 in Phillips et al. (198J}, we conclude that the 
pattern seen in Fig. [8^ is due to self-absorption (i.e., by a tem- 
perature gradient in the cloud core), rather than due to absorption 
by a colder (unrelated) foreground cloud. This is consistent with 
the fact that the CS-data reveal a higher temperature than the 
12 CO (Sect. 3.2.1). 

Again assuming LTE-conditions, we calculated the mass of 
the gas in the line wings from the 12 CO(2-l) data; the wings 
are not visible in the other two isotopomers. We adopted T ex - 
20 K, used T m \, (see Sect. 3.2.1), and considered both the opti- 
cally thick and -thin cases. The total velocity intervals, 24.7 - 
31.7 km s -1 (blue wing) and 36.7 - 43.7 km s (red wing), over 
which was integrated (in 1 km s~' -wide bins), were chosen such 
that the emission of the quiescent gas was excluded. 

In addition to the mass, we also calculated the flow's energy, 
momentum, size, and velocity (as outlined, e.g., by Lada 1985; 
Wouterloot & Brand [19991 ; the dynamical age was then derived 
from the ratio of the flow's size (~ 3.5 pc) and its velocity. 
The outflow parameters are collected in Table [2] Here we also 
give the flow's mechanical luminosity (L m = Energy/age) and the 
force needed to drive the flow (F m = Momentum/age). 

Our 12 CO(3-2) data, taken with a smaller beam size and 
smaller grid spacing, show self-absorbed profiles in the cloud 
core. Line wings are more extended than for CO(2-l), and reach 
20 kms -1 (blue) and 55 kms -1 (red), respectively. The top pan- 
els of Fig. [9] show a position-velocity diagram at A6-Q" (left), 
and a contour plot of the outflow lobes in the CO(3-2)-emission 
(right). 

We have convolved the CO(3-2) observations with a 
Gaussian beam to the resolution of the CO(2-l) data, and re- 
sampled all (2-1) and (3-2) spectra to a velocity resolution of 
0.55 km s" 1 . Using t] fss = 0.71, 0.80 for CO(3-2, 2-1), we have 
determined the line ratio R32 = r mb [CO(3 - 2)]/T mb [CO(2 - 1)] 
in each velocity channel where T m \, > 1.5<x for both lines. For 
optically thick lines the ratio R32 tends to 1 for large values of 
7kj n , while the limit for optically thin lines is 2.25 (Bachiller & 
Tafalla [1999l Codella et al. 120061 In the lower panels of Fig. [9] 
we show the results at two representative positions: near the cen- 
ter of the flow and the peak of the blue lobe (lower left) and at a 
more quiescent location in the cloud (lower right). Typical val- 
ues for R32 in the line center are 0.80-0.90, corresponding to 
7kin ~ 15 - 30 K, assuming both lines are optically thick; this 
is consistent with what we found from the CO(2-l) lines, and 



from a fit to the spectral energy distribution (see Sect. 3.2.3). In 
the line wings, typically R32 < 1.5, implying 7\i n < 40 K if 
both lines are optically thin away from the line center. But R32 
can reach values of 4 to 6 at some positions, as can be seen in 
the panel on the lower left in Fig. [9] This cannot be explained 
if both lines are optically thin or thick at the same time. On the 
other hand, the uncertainties in the high /?32-values are large, 
and assuming a more stringent criterion in their calculation, e.g., 
by requiring the value of T m \, to be at least 3<x in both transi- 
tions, would obviously remove most (though not all) points with 
anomalous values of R32. Ratios R32 of around 2 are found for 
optically thin lines and T^ a as 120 K, which is not unexpected 
in the shocked gas in molecular outflows (see, e.g., Hatchell et 
al. [1999 ; Hirano & Taniguchi 2001). A more accurate value of 
the temperature of the outflowing gas, and the possible existence 
of a gradient along the lobes, can be obtained only by observ- 
ing higher transitions and at higher resolution (the JCMT beam 
sizes at the frequencies of the CO(3-2) and (2-1) transitions cor- 
respond to 0.8 and 1.2 pc, respectively) 

Emission lines in the "tail" of the cloud, i.e., at offsets 
Aa - 60" (see Fig. [7J are less intense than in the core, with 
T mb < 5 K and < 7 K for CO(3-2) and (2-1), respectively. The 
lines are narrow (AVf w h m ~ 1.1-1.5 kms -1 ), becoming pro- 
gressively narrower towards the W-edge of the cloud. There is a 
slight velocity gradient in this part of the cloud, as can be seen in 
Figs. [8^ and[9](top, left). Line ratios R32 indicate typical values 
of Tkin ~ 8 + 2 K in this region. 

Following the same procedure as outlined above, we have 
calculated various outflow parameters from the CO(3-2) data; 
these too are collected in Table [2] The outflow in the CO(3-2) 
line is stronger and more extended in velocity, and more com- 
pact than that in the CO(2-l) line. The smaller size and higher 
average mass-weighted outflow velocity for the 3-2 flow result 
in a lower dynamical age. Other flow parameters are likewise in- 
fluenced by this, but as can be seen from Table [2] in general the 
results obtained from the two transitions are consistent within 
factors of a few. The lobes of this outflow have considerable 
overlap, indicating that the outflow axis is pointing close to the 
line-of-sight. As a consequence the size derived here is likely to 
be a lower limit, and age an upper limit. Values of F m and L m 
are thus lower limits too. 

In the calculation of the flow parameters in Table [2] we have used 
the local abundance of [CO]/[H2]. It is possible that the abun- 
dance ratio is lower at larger galactocentric distance. Gradients 
in C/H and O/H are similar (Wilson & Matteucci [19921 Rolleston 
et al. l2000l l. and assuming that the CO-abundance is determined 
by that in C, one expects [CO]/[H2] to be a factor of about 4 
lower in the far-outer Galaxy. This would increase all the param- 
eters listed in Table |2]by a factor of 4, except for the size, flow 
velocity, and age. Keeping this in mind, a comparison with out- 
flows associated with IR sources of similar luminosity (Beuther 
et al. 120021 Zhang et al. 2005) shows that this is not inconsis- 
tent with a typical outflow in an intermediate-mass star-forming 
region. 

In Fig. [9] top right panel, we show the CO(3-2) outflow 
superimposed on the K-image. Near offsets (+10", -3") and 
(-16",4"), bright extended spots are seen (cf. also Fig.[3h). The 
spots are brightest in K, and virtually absent in J (cf. Fig. QJ)). 
This emission might be due to the 2.1 fim line of H2, a manifes- 
tation of the interaction between outflow and ambient medium. 
This can be established only by taking a spectrum. 
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Fig. 8. a. Aff - velocity diagram at A6=0" of the 12 CO(2-l) emission. Contour levels are 0.25, 0.4(0.2)1(0.5)3(1)15 K. At the 
location of the IRAS point source (at 0,0) emission in the wings of the line profiles is evident, b. Map of the integrated wing 
emission. Blue (drawn): 24.7 < Vi sr < 3 1 .7 km s" 1 ; red (dashed): 36.7 < Vi sr < 43.7 km s . Contour levels and step are 0.5 Kkm s . 
The filled triangle marks the position of the peak of the SCUBA 850 /urn emission (— 1."9, 2."1); the observed positions are indicated 
with the crosses. 

Table 2. Parameters of the molecular outflow associated with WB89-789 1 



CO(2-l) 



Wing 




Mass 




Mom. 






Energy 




Size 


v out 


^dyn 




F m 






^111 








M G 




M G kms 


l 




10 44 erg 




pc 


kms~' 


10 5 yr 


10 4 


M kms 


yr 1 




10- 2 L G 




blue 


5.1 


3.2 


4.1 


18 11 


15 


6.7 


4.2 


5.4 


1.4 


3.5 


4.0 


0.5 


0.3 


0.4 


1.4 


0.9 


1.1 


red 


4.0 


2.6 


3.3 


14 9 


11 


4.7 


3.0 


3.8 


1.6 


3.3 


4.8 


0.3 


0.2 


0.2 


0.8 


0.5 


0.7 


CO(3-2) 


blue 


6.4 


4.0 


3.4 


28 18 


15 


15 


10 


8 


1.1 


4.4 


2.4 


1.2 


0.7 


0.6 


5.2 


3.3 


2.8 


red 


5.2 


3.3 


2.8 


27 17 


14 


20 


13 


11 


1.2 


5.1 


2.3 


1.2 


0.7 


0.6 


7.2 


4.6 


3.9 



t In case of three numbers per column, left: T cx = 20, t = 1, center: T cx = 20, r = 0, right: r ex = 40, r = 0. 
We used the local abundance ratio [CO]/[H 2 ] = 1.0 xl0~ 4 . 



Table 3. Cloud parameters derived from the continuum measure- 
ments. 





l i 


S 2 






pc 


iy 


M Q 


450 yum 


0.85 0.25 


12.9 




850 yum 


1.48 0.39 


1.8 




1.2 mm 


1.50 0.76 


0.5 




SED fit 
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1. r c , r\ p ; radius at 3<r-level and half-peak level, respectively; both 
corrected for beam size. 

2. Integrated flux density above the 3<x-level 

3. Assuming optically thin dust emission, a dust absorption coeff. 

1 cm 2 g-' at 250 GHz; the exponent of the freq. dependence of the dust 
optical depth /? = 1.8, and r,j ugt =23 K (see text for details). 



3.2.3. Dust continuum 

In Fig. [10] we have collected the maps of the dust continuum, 
made at 3 different wavelengths, and shown their location with 



respect to the core found in the molecular lines. The 13 CO emis- 
sion peaks very close to the maximum in the 850 fim map, while 
the C ls O peak is offset by about 10" (1 grid point ~ half a beam 
width) to the East. Flux densities at the three wavelengths were 
calculated by integrating all emission above the 3cr-contour, and 
sizes of the dust-emission cores were determined from the half- 
peak value contour and were corrected for beam size. The results 
are collected in Table [3] 

Figure QT| shows the spectral energy distribution (SED) of 
WB89-789. We have plotted our three FIR and mm-data points, 
in addition to the 4 IRAS flux densities, taken from the IPSC2. 
We have made a grey -body fit to the data points for 100 /mi 
onwards; the fit is shown as the drawn curve. The 12, 25, and 
60 yum points cannot be fit by the same curve. As shown by 
Molinari et al. d 19981 ) and Fontani et al. (120041) the flux densi- 
ties at these wavelengths usually arise from regions outside the 
core that are detected at the longer wavelengths, and one needs 
at least 2 grey-body fits to cover the entire SED. The SED- 
fit shown in Fig. [TJJwas obtained by assuming a dust opacity 
k v = k 25 o(v[GH z]/25Q'f, with k 25 o = 1 cm 2 g-' (Ossenkopf & 
Henning|l981). We obtained a good fit with fi — 1.8, and a dust 
temperature 7d us t = 23 K. The resulting dust mass is 11 M Q . 
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Fig. 9. Top. (Left) Aa - velocity diagram at A5=0" of the 12 CO(3-2) emission. Contour levels are 0.5, 1(1)15 K. (Right) Map of 
the integrated 12 CO(3-2) wing emission. Blue (drawn): 20 < Vi sr < 31.2 kms -1 ; red (dashed): 37.2 < Vi sr < 55 kms -1 . Contour 
levels and step are 2.0 Kkms -1 . The IRAS position is at (0,0). The white triangle with black outline marks the position of the peak 
of the SCUBA 850 //m emission. The underlying image is the K-frame of the region. Contours are drawn to highlight two extended 
features near offsets (+ 10", -3") and (- 16",4") that might mark spots where the outflow impacts with the ambient medium. Bottom. 
Line profiles and ratios at two positions: near the center of the flow and the peak of the blue lobe (left) and in a quiescent part of 
the cloud (right). The offsets are shown inside the panels. The drawn line profile is that of CO(3-2), the dotted one if for CO(2-l); 
the temperature-scale is shown on the right-hand axis. The open circles with error bars indicate the line ratio R32 for those velocity 
channels where T m ^ > \.5cr in both lines; values are shown on the left-hand axis. Short-dashed and long-dashed lines mark the 
upper limit of the line ratio if both lines are optically thick or thin, respectively (see text). 



For a gas-to-dust ratio of 100, the total mass of the core detected 
in the (sub-)mm continuum is 1 100 M Q , consistent with what is 
found from the 13 CO and C 18 0-data (see Table [B- We should 
keep in mind, however, that the dust properties and the gas-to- 
dust ratio in the FOG may differ from those found near the Sun. 

4. The distance to WB89-789 

The presence of an outflow, together with the H2O maser emis- 
sion detected by Wouterloot et al. ( 1993) demonstrates that star 
formation is still active in this object. If the kinematic distance of 
11.9 kpc can be confirmed, this would be the star-forming region 
furthest from the Galactic center yet found in the Galaxy. The 
distance was derived from the 12 CO velocity, using the Brand & 
Blitz ( 119931 1 rotation curve. WB89-789 is at a galactic longitude 
of 195°82; the determination of kinematic distances is notori- 
ously uncertain in directions near the (anti-) center, where Vi sl is 
not very sensitive to distance, and random motions become rel- 
atively influential in the distance determination. If we assume 
a cloud-cloud velocity dispersion of ~5 kms -1 (see Brand & 
Blitz 1993), then the range of galactocentric distance of WB89- 
789 is 17 < R < 25 kpc (with a corresponding range in heliocen- 
tric distance 8.6 Hd'H 17.1 kpc). 



As shown in Fig. [5^, star 33 (indicated in Fig. |3J, if it is a 
main sequence star, would be its earliest type member that does 
not have NIR excess, and would have a spectral type ~B0, if the 
cluster is at a distance of 11.9 kpc. To check this, we have ob- 
tained a spectrum of this star (see Fig.fTZl). From the way the con- 
tinuum rises towards longer wavelengths it is clear that star 33 
is not an early-type star. We have compared the TNG-spectrum 
with those in the library of Jacoby et al. ( I1984I ). To reduce the 
noise we smoothed the spectrum to a resolution of 5.4A, com- 
parable to the 4.5A-resolution of the Jacoby et al. ( 1984) spec- 
tra. From a visual inspection we conclude that the spectrum of 
star 33 agrees best with that of a K3 Ill-star (Jacoby et al. ( 1984) 
spectrum No. 100; this spectrum is also shown in Fig.[T2l. This 
implies that star 33 has a visual extinction Ay ~ 5.5, and a geo- 
centric distance of 10.7 kpc, i.e., R =s 19.0 kpc. Identification 
as a K3 Ill-star is consistent with the location of star 33 in the 
color-magnitude diagram of Fig. [5] dereddening brings it first on 
the giant branch, at around spectral type K3. 

This result is rather surprising: star 33 is embedded in neb- 
ulosity (cf. Fig. [3^), making it a likely member of the cluster 
associated with the IRAS source. On the other hand, a K3III- 
star would originate from a 9-10 M ms-star (i.e., spectral type 
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Fig. 10. Images of continuum emission at 450, 850, and 1200 //m, and overlays of the 850 /mi emission with the 13 CO and C ls O 
emission. The IRAS position is at (0,0). a. SCUBA 450 //m map. Contour values 1(1)5 Jy/beam. The white contour indicates the 
half-peak level, b. SCUBA 850 //m map. Contours 0.1(0.1)0.7 Jy/beam. White contour as in a. c. SIMBA 1.2 mm map. Contours 
45(30)230 mJy/beam. White contour as in a. d. SCUBA 850 /mi map with C ls O(2-l) contours superimposed. Contour values are 
0.3(0.3)2.1 Kkms -1 . The white contour indicates the 850 //m emission half-peak level, while the triangle marks its peak. e. As d, 
but with 13 CO(2-l) contours superimposed. Contour values 1(2)15 Kkms -1 . White contour and triangle as in d. 
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Fig. 11. Spectral energy distribution (SED) of WB89-789. 
Along the axes are \ogA and logf\ The data shown are our 
own measurements of the continuum emission at 450, 850, and 
1200 /mi, together with the flux densities in the 4 IRAS bands 
(IPSC2). The curve shows the best fit to the data points at wave- 
lengths > 100 //m, which has 7d us t = 23 K, and a dust absorption 
coefficient of 1.8, and which was obtained by assuming a dust 
opacity at 250 GHz of 1 cm 2 g _1 . 



B2-B3), and has an age of about 20-25 Myr. It is possible that 
star formation has been going on in this region for a much longer 
time, and that the IRAS source, the water maser, the outflow, and 
the NIR-excess stars are just the latest episode. The earlier oc- 
currence of star formation (giving rise to star 33) would have 
removed much of the original molecular cloud, perhaps explain- 
ing why the cloud in which star formation is currently taking 
place is so small. 

If this scenario is correct, then it still has to be explained how 
star formation can be sustained over such a long time, especially 
in this part of the Galaxy. As mentioned in the Introduction, if 
star formation is to be initiated by external triggers, then it is dif- 
ficult to see how this can work efficiently in the far-outer Galaxy, 
where external triggers (SN explosions, spiral arms) are much 
more scarce or weak. On the other hand, as already found by 
Wouterloot et al. (1988) and Snell et al. (2002), once star for- 
mation does occur in FOG-clouds, it appears to proceed as it 
does locally. It should also be pointed out that there are other 
examples of older stars located in younger clusters. Kumar et al. 
(2004) found three late-type giants in two young (1.5-1.9 Myr) 
open clusters, indicating that star formation has been going on 
there for about 25 Myr. However, the spectral classification of 
star 33 is very tentative, and further speculation should wait un- 
til a better spectrum is available. 



5. Conclusions 

We have studied the star-forming region associated with 
IRAS06 145 +1455 (WB89-789) in detail. Its location near the 
edge of the molecular disk (at a galactocentric distance of 
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Fig. 12. Spectrum of star 33 (see Fig.0 taken with DOLORES 
at the TNG. The spectrum has been smoothed to a resolution 
of ~ 5 .4 A. Also shown is the spectrum of a K3 Ill-star (from 
Jacoby et al. [T984 arbitrarily shifted vertically to allow com- 
parison). Because of the poor wavelength calibration, the spec- 
trum of star 33 has been shifted in A to get correspondence be- 
tween the features in both stars. The amount of shifting neces- 
sary agrees with that which we determined from a comparison 
of the lines in our flux calibrator (Hiltner 600) and Jacoby et al. 
( 1984) star No. 13, which has a similar spectral type. 

20.2 kpc) makes this the galactic star forming region furthest 
from the Galactic center found up to now. An embedded cluster 
containing about 60 stars was found in JHK-images. The radius 
of the cluster is about 1.3 pc, within which the average stellar 
surface density is 12 pc~ 2 ; within a core radius of 0.4 pc the 
stellar surface density is 33 pc~ 2 . The stars have a foreground 
extinction of about 6 mags and an internal extinction of up to 
30 mags, though typically ~ 7.5 mags. At least 14 of the em- 
bedded stars show NIR-color excess in a J-H vs. H-K diagram, 
suggestive of the presence of disks around these stars; three of 
these are possibly of Class I, while the locations of the others 
suggest they are Class II objects. 

An optical spectrum of what seemed like it could be the earli- 
est type cluster member with normal extinction shows a tentative 
likeness to that of a K3 III star. The distance of this star would 
then be 10.7 kpc (galactocentric distance 19.0 kpc), consistent 
with the kinematical distance of the region. The presence of this 
star in the cluster suggests that star formation has proceeded for 
the past several 10 7 years. The absence of radio continuum emis- 
sion towards this source, at the sensitivity of our observations, 
indicates that no star of type B0. 5 V or earlier is present. 

JCMT 13 CO and C lg O and IRAM CS observations show 
that the cluster is embedded in a ~ 500 - 1500 M Q (M vir ) 
core. A JCMT 12 CO map shows that the core is contained in 
a ~ 5000 M cloud, consisting of two clumps. The clump with 
the embedded cluster shows the existence of moderate outflow 
emission. (Sub-)mm continuum maps (JCMT and SEST) at 0.45, 
0.85, and 1.2 mm show the presence of a dust core centered at 
the cluster. With standard assumptions the dust mass is 10 M ; 
assuming the canonical gas/dust ratio of about 100, the total core 
mass is ~ 1000 M , which is consistent with what is found from 
13 CO and C ls O data. The question now is which mechanism 
started the star formation process in this cloud, and how it can 
be sustained over such a long time (i.e., at least 20 - 25 Myr, the 
age of a K3 Ill-star). 
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Table 4. NIR photometry data 
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Table 4. continued 



nr 


a(2000) 
hm s 


o ' " 


K 






0" H 


j 


V J 


H-K 


J-H 


46 


06: 


17:23.79 


+ 14:54:27.9 


16.25 


0.10 


16.74 


0.08 


17.78 


0.12 


0.487 


1.037 


47 


06: 


17:23.64 


+ 14:54:36.3 


15.74 


0.09 


16.26 


0.07 


17.45 


0.12 


0.528 


1.190 


48 


06: 


17:23.56 


+14:55:15.8 


17.32 


0.16 


17.70 


0.12 


18.68 


0.23 


0.378 


0.978 


49 


06: 


17:23.56 


+14:54:52.8 


17.06 


0.16 


18.32 


0.21 


19.60 


0.52 


1.269 


1.280 


50 


06: 


17:23.54 


+ 14:54:50.1 


16.27 


0.08 


16.96 


0.09 


18.12 


0.15 


0.691 


1.157 


51 


06: 


17:23.53 


+14:54:39.6 


16.15 


0.08 


16.50 


0.09 


17.71 


0.11 


0.348 


1.212 


53 


06: 


17:23.46 


+14:55:06.6 


16.96 


0.13 


17.39 


0.12 


18.83 


0.29 


0.433 


1.433 


54 


06: 


17:23.42 


+14:54:40.0 


16.19 


0.07 


18.29 


0.26 


19.34 


0.43 


2.097 


1.053 


56 


06: 


17:23.30 


+ 14:55:14.9 


15.42 


0.05 


15.60 


0.11 


15.95 


0.14 


0.173 


0.350 


57 


06: 


17:23.32 


+ 14:54:23.0 


16.14 


0.09 


16.75 


0.06 


17.80 


0.11 


0.605 


1.057 


60 


06: 


17:23.23 


+14:54:28.2 


15.52 


0.08 


16.25 


0.07 


17.57 


0.13 


0.730 


1.318 


62 


06: 


17:23.14 


+ 14:54:36.9 


16.25 


0.08 


16.79 


0.06 


18.02 


0.13 


0.544 


1.224 


63 


06: 


17:23.08 


+ 14:54:31.8 


17.54 


0.24 


18.86 


0.42 


19.63 


0.54 


1.325 


0.762 


64 


06: 


17:23.02 


+14:54:58.0 


13.68 


0.07 


13.97 


0.07 


14.58 


0.22 


0.288 


0.614 


65 


06: 


17:23.05 


+ 14:54:22.0 


16.52 


0.11 


17.01 


0.06 


18.14 


0.14 


0.487 


1.130 


66 


06: 


17:22.82 


+ 14:55:13.1 


15.74 


0.06 


16.33 


0.05 


17.33 


0.07 


0.581 


1.006 


68 


06: 


17:22.75 


+ 14:54:59.5 


14.12 


0.06 


14.59 


0.12 


15.70 


0.08 


0.472 


1.110 


69 


06: 


17:22.70 


+ 14:54:41.6 


13.09 


0.06 


13.63 


0.04 


14.67 


0.03 


0.542 


1.045 


70 


06: 


17:22.64 


+14:54:56.4 


14.99 


0.06 


15.19 


0.03 


15.88 


0.13 


0.201 


0.693 


71 


06: 


17:22.64 


+ 14:54:33.8 


15.72 


0.05 


15.89 


0.09 


16.40 


0.14 


0.166 


0.512 


72 


06: 


17:22.45 


+ 14:54:55.0 


15.79 


0.07 


16.23 


0.05 


17.28 


0.08 


0.435 


1.052 


76 


06: 


17:25.10 


+ 14:54:22.4 


15.89 


0.08 


16.21 


0.09 


16.78 


0.07 


0.327 


0.562 


80 


06: 


17:24.25 


+ 14:54:44.2 


16.58 


0.23 


17.18 


0.13 


18.10 


0.15 


0.607 


0.915 


81 


06: 


17:24.28 


+14:54:26.6 


16.76 


0.25 


17.79 


0.25 


18.38 


0.25 


1.027 


0.593 


83 


06: 


17:23.95 


+ 14:54:41.8 


14.74 


0.28 


17.02 


0.27 


18.11 


0.19 


2.287 


1.091 


86 


06: 


17:23.68 


+14:55:18.3 


17.03 


0.22 


17.67 


0.21 


18.94 


0.40 


0.639 


1.275 


87 


06: 


17:23.37 


+14:54:42.2 


17.75 


0.32 


17.83 


0.16 


19.22 


0.40 


0.085 


1.384 


88 


06: 


17:23.13 


+14:54:26.6 


16.24 


0.19 


16.49 


0.08 


17.13 


0.09 


0.254 


0.638 


89 


06: 


17:23.09 


+14:54:57.2 
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17.28 


0.16 


18.81 


0.47 


0.739 


1.535 
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15.11 


0.07 


15.22 
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0.04 
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1.220 


1.913 
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0.08 


16.35 


0.06 
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0.18 


17.12 


0.23 


17.84 


0.21 


0.541 
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17:23.73 
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0.25 


17.46 


0.32 
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0.18 
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